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The implementation of the SSG/LRR-u; differential Reynolds stress model into the 
NASA flow solvers CFL3D and FUN3D and the DLR flow solver TAU is verified by study- 
ing the grid convergence of the solution of three different test cases from the Turbulence 
Modeling Resource Website. The model’s predictive capabilities are assessed based on four 
basic and four extended validation cases also provided on this website, involving attached 
and separated boundary layer flows, effects of streamline curvature and secondary flow. 
Simulation results are compared against experimental data and predictions by the eddy- 
viscosity models of Spalart-Allmaras (SA) and Menter’s Shear Stress Transport (SST). 
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Cartesian component of diffusion tensor, m 1 2 /s 3 
Coefficient of generalized gradient diffusion model 
Coefficient of simple gradient diffusion model 
Distance to closest wall, m 
Menter’s blending function 
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h Equivalent grid spacing 

k Specific kinetic turbulence energy, m 2 /s 2 

H Height of jet orifice or bump or square duct, m 

L Unit length or flat plate length, m 

M Mach number 

N Grid point number 

Pij Cartesian component of Reynolds stress production tensor, m 2 /s 3 

Rij Cartesian component of specific Reynolds stress tensor, m 2 /s 2 

R~ Non-dimensional Cartesian component of specific Reynolds stress tensor in inner scaling 

R xy Specific Reynolds shear stress 

Rei Reynolds number based on length £, with £ = c, D, H , or L 

Sij Cartesian component of strain rate tensor, 1/s 

S*j Cartesian component of traceless strain rate tensor, 1/s 

T Temperature, K 

t Time, s 

U Streamwise Cartesian velocity component, m/s 

Uq Streamwise Cartesian velocity component at jet center, m/s 

U\ Streamwise Cartesian velocity component, m/s 

U i Cartesian component of velocity vector, m/s 

Uoo Free stream velocity, m/s 

u T Friction velocity, m/s 

Wij Cartesian component of rotation tensor, 1/s 

w{y/8) Wake-function 

x Streamwise Cartesian co-ordinate, m 

Xi Cartesian co-ordinate, m 

y Wall-normal Cartesian co-ordinate, nr 

y' Wall-normal Cartesian co-ordinate with origin at the surface, m 

y + Non-dimensional wall-normal Cartesian co-ordinate in inner scaling 
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Subscripts 

Ref Reference conditions 
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Incidence, deg. 

Coefficient of ui production 
Coefficient of ui destruction 
Boundary layer thickness, m 
Displacement thickness, nr 
Kronecker symbol 
Isotropic dissipation rate, nr 2 /s 3 

Cartesian component of dissipation rate tensor, nr 2 /s 3 
Argument to Menter’s blending function F\ 

Von Karman constant 

Molecular (dynamic) viscosity, kg/ (m • s) 

Molecular (kinematic) viscosity, m 2 /s 2 
Wake parameter 

Cartesian component of pressure-strain correlation, m 2 /s 3 
Density, kg/m 3 

Coefficient of to cross-diffusion term 
Coefficient of lo diffusion term 
Wall shear stress, kg/(m-s 2 ) 

Specific dissipation rate, 1/s 
Simple averaged value 
Mass weighted averaged value 
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I. Introduction 


Computational Fluid Dynamics (CFD) is routinely applied in contemporary aircraft design. Nevertheless, 
the predicted result of any flow simulation depends on the accuracy of the underlying representation of the 
flow physics, i.e. , the effect of turbulence, and the accuracy of the numerical method solving the corresponding 
equations. 

While the numerical error, at least in principle, can be reduced by refining the computational mesh 
and, in time-accurate simulations, the time step, reducing the modeling error is much more demanding. 
Theoretically, a turbulence model can be completely avoided by performing a Direction Numerical Simulation 
(DNS), which then requires resolving all turbulent fluctuations down to the smallest scales. Alternatively, in 
a Large Eddy Simulation (LES), the resolution requirements can be relaxed, but at the expense of modeling 
the effect of the smallest scales. Unfortunately, the size of the smallest scales to be resolved decreases rapidly 
with increasing Reynolds number, currently restricting the application of such scale resolving methods. 1 

For this reason, turbulence models based on the Reynolds-averaged Navier-Stokes (RANS) equations are 
still the backbone of industrially applied CFD methods (for an overview see Ref. 2). Within this framework, 
only the mean flow is considered, where the net effect of turbulence enters into the momentum balance via 
additional Reynolds stresses. As a consequence, these additional unknowns need to be expressed in terms of 
known mean quantities, which constitutes the art of RANS turbulence modeling. 

The most widespread RANS modeling approach is based on the Boussinesq hypothesis, assuming the 
turbulence to effectively increase the viscosity of the fluid by a flow dependent eddy viscosity. 2 Thus, only 
one single parameter needs to be provided that can be calibrated to correctly predict (for example) the 
turbulent shear stress in a boundary layer. Indeed, eddy-viscosity models (EVM) like the Spalart-Allmaras 
(SA) model 3 or the Menter Shear Stress Transport (SST) model 4 have proven to give reliable predictions in 
a wide range of applications and are therefore heavily used within the aeronautics industry. 

However, pushing the applications to the borders of flight, the accuracy of the applied eddy-viscosity 
models degrades, where the prediction of flow separation appears to be the greatest challenge. One possibility 
for improvement consists in replacing the Boussinesq hypothesis by solving the individual transport equations 
for the Reynolds stresses. Corresponding Second-Moment Closures (SMC), also known as Reynolds stress 
models (RSM), allow for a more detailed representation of the flow physics, but nevertheless require modeling 
of individual terms in the underlying transport equations (for an overview see Ref. 5). Hence, as with eddy- 
viscosity models, the accuracy of Reynolds stress models needs to be assessed individually. 

The SSG/LRR-w Reynolds stress model 6 has been demonstrated to be applicable to fairly complex 
aeronautical flow problems. It has therefore been selected by NASA and DLR for a common verification and 
validation exercise based on the test cases on the Turbulence Modeling Resource (TMR) Website. 7 


II. SSG/LRR-w Model 


The SSG/LRR-w model has been developed within the EU-project FLOMANIA 8 and, as any Reynolds 
stress model, is based on the Reynolds stress transport equation 
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where p is the mean density, U t represents the Cartesian components of the mean velocity vector and ~pRij 
the Cartesian components of the Reynolds stress tensor. 

The Cartesian components of the production term are given by 
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and do not need modeling, since all quantities are provided by the system of equations to be solved. 

The Cartesian components of the pressure-strain correlation are formally modeled according to Speziale, 
Sarkar and Gatski (SSG) 9 as 
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where k = Ru/2 is the specific kinetic turbulence energy and 
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denote the Cartesian components of the simple strain rate tensor, the traceless strain rate tensor and the 
rotation tensor, respectively. The C, and C* are model dependent coefficients. 

Dissipation is modeled as an isotropic tensor with Cartesian components 
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(8) 


where the isotropic dissipation rate e is provided indirectly by a transport equation for the specific dissipation 
rate to given further below. 

Two different models are available representing non-convective fluxes, simple gradient diffusion 10 
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where p is the mean dynamic fluid viscosity and D^ sd ' ) and D represent respective model coefficients. 
Finally, Menter’s baseline equation 4 
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is employed for providing the isotropic dissipation rate 

e = C^kx, (12) 

where C ^ = 0.09. As with the SST model, the coefficients (f> = a u , /3 W , , oq are blended according to 

<j> = + (1 - FJ <j>^ (13) 

between the bounding values associated with the w-equation by Wilcox 12 near walls (Fi = 1) and the 
standard e-equation 13 at the boundary layer edge (iq = 0), where the blending function is given by 4 
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Table 1. Bounding Values of ^-Equation Coefficients 


Q-oj fiuj @d 

<^ £ ) 0.44 0.0828 0.856 1.712 

0.5556 0.075 0.5 0 


based on the wall-distance d. The corresponding values of the coefficients are given in Table 1. 

The fundamental idea of the SSG/LRR-w model is to combine the e-based SSG model for the pressure- 
strain correlation with the simpler model by Launder, Reece and Rodi (LRR) 14 near walls in the w-based 
formulation of Wilcox. 15 Since the LRR-model can be written as a subset of the above SSG-model, Eq. (3), 
the same blending of the Reynolds stress model coefficients C*, C* , D^ sd ^ and D < ' GD ' ) can be applied as to the 
coefficients of the w-equation. Thus, all coefficients change consistently from LRR+w near walls to SSG+e 
at the boundary layer edge. The corresponding bounding values are given in Table 2. Note that the value of 
the coefficient c£ RIR> as well as the generalized diffusion coefficient D ( GD ^ have been re-calibrated for better 
agreement with the log-law velocity profile in a zero pressure gradient boundary layer. 16 Therefore, on the 
Turbulence Modeling Resource Website, 7 the model is denoted SSG/LRR-RSM-w2012 when employing the 
generalized diffusion model and SSG/LRR-RSM-w2012-SD when employing the simple gradient diffusion 
model. In this paper only the generalized gradient diffusion model is used so that for brevity the model is 
termed SSG/LRR-w only. 
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III. Verification of Implementation 


III. A. Method 

In this article, verification refers to the objective that a turbulence model, here the SSG/LRR-w, is im- 
plemented correctly, i.e., according to the specified model equations, including boundary conditions. 7 Since 
giving a strict proof might be rather difficult, a more heuristic weak form of verification is applied, comparing 
simulation results obtained with independently developed flow solvers on a sequence of successively refined 
grids. Good agreement between different solutions in the limit of grid independence enhances the confidence 
in the correct implementation of the model into all of the tested codes. As already said, the comparison will 
be restricted to the 2012 version of the SSG/LRR-w model with generalized gradient diffusion. 

III.B. Flow Solvers 

The SSG/LRR-w model has been implemented independently into three different production codes of NASA 
and DLR: CFL3D, FUN3D and TAU, solving the compressible Reynolds-averaged Navier-Stokes equations 
together with the respective turbulence equations. CFL3D is a cell-centered structured-grid code developed 
at NASA, whereas FUN3D and TAU are node-centered unstructured-grid codes developed by NASA and 
DLR, respectively. All three are based on a finite- volume discretization of formally second order accuracy in 
space, where the preferred specific discretization scheme varies between the different codes and organizations. 
The advection terms of the turbulence equations are usually solved with first order accuracy only. Most cases 
to be shown below used this option. However, for a few of the cases, TAU also employed second order accuracy 
for the turbulence advection terms to explore their influence. 
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III.C. Results 


Three two-dimensional verification cases have been selected from the TMR Website:' the flow over a flat 
plate at zero pressure gradient, the plane jet, and the flow over a bump. These cases are described in detail 
in the following sections. 

III.C.l. Flat plate at zero pressure gradient 

The flow over a flat plate is computed with the mean flow boundary conditions defined according to the 
TMR Website, Fig. I. 7 The Mach number is M = 0.2, Reynolds number is ReL = 10 • 10 6 based on the 
flat plate length L 1 and reference temperature is Tfl e / = 300 K. The computations have been carried out 
on the sequence of five structured grids, also provided on the TMR Website, 7 ranging from TV = 35 x 25 
points on the coarsest level to TV = 545 x 385 points on the finest level, where the coarser levels have been 
generated by omitting every second grid line in each coordinate direction. The TAU computations have 
been performed with two different spatial discretization schemes for the advection terms of the turbulence 
equations, allowing us to study the influence of this term on the overall numerical accuracy. 

Figure 2 (left) shows the variation of the predicted total drag coefficient Cd with grid resolution in terms 
of an equivalent grid spacing h = 1/ \/N. The differences between the predicted results reduce significantly, 
when refining the mesh, indicating one common solution in the limit of infinitely small cells. This is further 
verified by Fig. 2 (right), depicting the dependence of the local skin friction coefficient, c/, at x/L = 0.485 
with the grid resolution. The differences between the different predictions on the finest mesh appear even 
smaller than for the total drag coefficient. Note that in agreement with boundary layer theory, 17 in this 
particular case, no advantage of second order accuracy for the advection terms can be observed on the finer 
grid levels. 

As it turns out, the somewhat larger differences in total drag are due to differences in the vicinity of 
the flat plate leading edge only. This is revealed by Fig. 3, showing a log-log plot of the local skin friction 
coefficient over the flat plate predicted on the finest mesh. Larger differences are only observed up to 
x/L Ki 5 • 10 -3 , which can be attributed to the singularity at the leading edge. 

The very good agreement between the different predictions on the finest mesh is confirmed, when com- 
paring profiles of various quantities normal to the wall. Figure 4 shows the predicted velocity profiles at 
x/L = 0.485 normal to the wall in inner scaling. The inner scaling is defined by relating all quantities to the 
kinematic viscosity, u, of the fluid and the friction velocity 

(16) 

where t w is the local shear stress at the wall and p is the fluid density. Corresponding quantities are indicated 
with a superscript + , i.e. , y + = yu T /v , u + = U\/u T , R/j = Rij/v%, where the tilde symbol over U\ and Rij 
has been dropped for convenience. Virtually no difference can be observed between the predicted velocity 
profiles over the complete boundary layer height. 

Figure 5 shows the corresponding profiles of the non-dimensional Reynolds stresses in inner scaling at 
the same position, x/L = 0.485. Small deviations can be observed only towards the boundary layer edge; 
in particular, the second-order turbulence advection terms steepen the gradients in this region. However, 
the agreement between the CFL3D and TAU results (both with first-order turbulence advection) is nearly 
perfect. 

III.C. 2. Plane jet 

In this case, a plane jet is developing behind an orifice, from where fluid is entering into an ambient stream 
at lower uniform velocity. Corresponding computations have been carried out on the series of five structured 
meshes provided on the TMR Website 7 using the boundary conditions as depicted in Fig. 6. Reference Mach 
number is M^ef = 0.5, Reynolds number is Ren = 50, 000 based on height of the orifice, and reference 
temperature is Tn e f = 300 K. Note that the orifice is modeled as an infinitely thin plate subjected to the 
flow on both sides. 

The original grids are composed of three blocks with a total point number of TV = 1, 280 on the coarsest 
up to TV = 327, 680 on the finest level. The coarser grids are obtained from the respective next finer one by 
omitting every second grid line in both directions. 
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Figure 7 (left) shows the variation of the total drag coefficient Cd on the thin plate with the grid resolution 
in terms of an equivalent grid spacing h = 1/y/N. As to be expected, the differences between the predicted 
values decreases with increasing mesh resolution. Nevertheless, even on the finest mesh there remains a very 
small difference. 

Additionally, Fig. 7 (right) presents the non-dimensional streamwise velocity, Uo/aaef (where a# e / is the 
reference speed of sound), on the centerline of the jet versus the equivalent grid spacing h at three different 
positions x/H downstream of the orifice. Once again, the differences between the predictions reduce with 
increasing mesh resolution, but the variation differs between the three positions. Interestingly, the largest 
differences remain at x/H = 29.2468, which is located between the two other locations. 

That the remaining differences are indeed very small can be seen from Fig. 8 showing the development 
of the non-dimensional streamwise velocity on the jet centerline along the flow computed on the finest mesh. 
On this scale, there is virtually no difference between the predictions by the different codes. 

This observation is confirmed when comparing profiles of the non-dimensional velocity and Reynolds 
stresses across the jet. Figure 9 shows the corresponding results for the non-dimensional velocity profile on 
the finest mesh at position x/H = 29.2468, where the largest differences have been observed in the centerline 
velocity. All results virtually collapse on a single line. A close-up of the velocity near the centerline reveals 
small differences that may reflect differences in the numerical treatment of symmetry planes in the different 
flow solvers. Similarly the non-dimensional Reynolds stress profiles, shown in Fig. 10, nearly agree perfectly. 

III.C.3. Two-dimensional bump in a channel 

The two-dimensional flow over a bump in a channel is a verification test case on the TMR Website 7 aimed 
at an increased complexity of the flow field due to the variation of the pressure gradient over the bump. The 
contour of the lower channel wall is defined by 7 

! 0 x/H < 6 

{«n[f (^f -l)]} 4 6 < x/H < 24 , (17) 

0 x/H > 24 

where H is the height of the bump. It is shown in Fig. 11 together with the corresponding boundary 
conditions. The Mach number is M = 0.2, Reynolds number Ren = 150,000 based on the bump height, 
and reference temperature is = 300 K. 

Computations have been carried out on the series of five structured meshes provided on the TMR Website 7 
with point numbers ranging from N = 89 x 41 on the coarsest to N = 1,409 x 641 on the finest level. As 
usual, the coarser grids have been obtained by omitting every second grid line of the respective next finer 
level. 

Figure 12 shows the variation of different total force coefficients with the equivalent grid spacing h = 
1/y/N. In particular, the differences in the lift coefficient Cl, Fig. 12 (left), and the pressure drag coefficient 
Cd, p , Fig. 12 (right), are decreasing with increased mesh resolution. For the viscous drag coefficient Cd,v, 
Fig. 12 (right), the effect is less prominent, in particular when plotted to the same scales as the pressure 
drag coefficient. Nevertheless, Fig. 13 shows that the local skin friction coefficient at three positions around 
the crest of the bump converges to almost identical values with all three codes. 

Figure 14 shows the distributions of the local pressure coefficient C p (left) and the local skin friction 
coefficient c/ (right) along the viscous wall, computed on the finest mesh with N = 1, 409 x 641 points. The 
results are virtually identical, except for some differences in c/ near the inflow plane that vanish rapidly. 

Figure 15 shows the profile of the streamwise velocity in inner scaling at x/H = 24, corresponding to 
the end of the hump. At this position the complete history effects of the flow over the bump are expected 
to be present. Compared to the flow over a flat plate at zero pressure gradient, Fig. 4, the velocity profile 
is largely distorted, and the corresponding Reynolds stresses in Fig. 16 show a significant peak close to the 
boundary layer edge. Nevertheless, all predictions from the three codes virtually collapse to a single result 
except for some insignificant differences near the Reynolds stress extrema. 

III.D. Summary of Verification 

In all three considered verification cases, flat plate, plane jet and bump flow, the predictions of CFL3D, 
FUN3D and TAU apparently converge towards the same solution. While integral indicators, like total force 
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coefficients, still show some very small remaining differences on the respective finest mesh, there is hardly 
any difference visible when plotting skin friction or pressure distributions in standard scales. Furthermore, 
very good agreement has been observed in velocity and Reynolds stress profiles. 

Therefore, there is high confidence that the SSG/LRR-tu model has been implemented into all three 
codes correctly according to its specification. Under the assumption of a sufficiently fine grid and identical 
boundary conditions, predictions with any of these codes should yield the same results for this turbulence 
model. On arbitrary grids, however, differences in the results may be observed due to differences in the 
numerical methods. For this reason, results for the following validation cases will be cross-checked as far as 
possible. 


IV. Validation 

In this article, validation refers to the question of how well real flows are represented by a turbulence 
model, e.g., the SSG/LRR-w model. Therefore, in this section, the focus is on comparisons of numerical 
simulation results with experimental data. According to the TMR Website, 7 the different test cases will 
be separated into basic and extended ones, where the focus will be on wall-bounded flows. Based on the 
experience presented on the TMR Website, computations on the second finest grid level are usually considered 
sufficient. SSG/LRR-w model predictions with TAU, CFL3D, and/or FUN3D will be presented to get an 
indication of any possible remaining numerical influence. Additionally, results obtained with CFL3D using 
the well-established Spalart-Allmaras (SA) model 3 and the Shear Stress Transport (SST) model by Menter 4 
are included for comparison. The verification of these models in CFL3D is documented on the TMR Website. 

IV. A. Basic Validation Cases 

IV.A.l. Flat plate at zero pressure gradient 

The flow over a flat plate at zero pressure gradient has also been selected as a basic validation case, where the 
fine grid results of the above verification study are used. The flow conditions refer to a Reynolds number of 
ReL = 10 • 10 6 based on the total length L of the flat plate, corresponding to the experimental conditions of 
Wieghardt. 18 According to the TMR Website, the Mach number is set to M = 0.2 and reference temperature 
is 300 K. 7 Numerical results for the SSG/LRR-cu model on the fine grid are taken from the above verification 
study and are complemented by fine grid results obtained with CFL3D, using the SA and the SST models. 

Figure 17 (left) shows the predicted velocity profiles in inner scaling, extracted at x/L = 0.95, compared 
to the analytical expressions for the viscous sublayer 

u + = y + , (18) 


valid up to y + w 5, and the so-called log-law 

u + = —\ny + + C (19) 

K 

with standard values for the von Karman constant n = 0.41 and the intercept C = 5. The importance of 
the log-law lies in the fact that it is apparently valid in any wall-bounded incompressible turbulent flow, 
independent of the Reynolds number. 19 There is close agreement between the SST and the SSG/LRR- 
l o predictions over almost the whole boundary layer height. The SA model predicts a smaller buffer layer 
between the viscous sublayer and the log-law region, which is probably in better agreement with experimental 
observation, although not shown here. Closer inspection of the region 100 < y + < 1000, depicted in Fig. 17 
(right), reveals that all models predict a slightly steeper velocity gradient than would follow from the log-law 
with k = 0.41. 

Figure 18 shows the skin friction distribution along the flat plate predicted by the different turbulence 
models compared to the experimental results by Wieghardt 18 as provided by Coles and Hirst. 20 All predic- 
tions are close to one another and, at least for lower local Reynolds number, tend to slightly underestimate 
the Wieghardt data. Closer inspection of the high Reynolds number region, Fig. 18 (right), reveals that 
the SSG/LRR-tu result virtually coincides with the SST result, whereas the SA model gives slightly higher 
skin friction, which apparently is in better agreement with the Wieghardt data. Nevertheless, the differences 
seem to be small, and all three turbulence model results appear acceptable. 
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Note that the skin friction is an indicator for the boundary layer thickness, 5, and correspondingly the 
displacement thickness, S * , as can be seen from the analytical velocity profile suggested by Coles 21 


u + = — In y + + C H w(y/S ), 

K K 


( 20 ) 


where w(y/5) is the wake function with the limiting values 

(21) 
( 22 ) 


w(y/8 = 0) = 0, 

w{y/8 = 1) = 2 


and the integral 


w d(y/8) = 1. 


(23) 


Jo 

The wake function is supposed to be universal, but is scaled with the flow dependent wake parameter II, 
assuming a value of II = 0.55 for the boundary layer of a flat plate. 21 The Coles profile, Eq. (20), describes 
the velocity profile in the whole boundary layer, except in the viscous sublayer and the buffer layer near 
walls, i.e., below y + « 30-50. 

Evaluating the velocity profile, Eq. (20), at the boundary layer edge, y = S, and noting that = 
Uoo/u T = \J 2,/cf (where oo represents the farfield reference conditions), one can re-write Eq. (20) to obtain 
an equation for the boundary layer thickness non-dimensionalized with the flat plate length L 


S_ 
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Re L 




{kC + 211) } , 


(24) 


where ReL is the Reynolds number based on the flat plate length. Clearly, the larger the skin friction 
coefficient, c/, the smaller the boundary layer thickness, S. As shown by Coles, 20 the displacement thickness, 
6*, can be related to the boundary layer thickness, 8 , by 


8 * 1/00 

Su T 



d{y/S) 


1 + 11 
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(25) 


ignoring the deviation of the velocity profile, Eq. (20), from the viscous sublayer relation near the wall. After 
some re-arrangement one obtains 


6 * 1 1 + n / 2 

T'ife — exp {' c V^-( n) J ’ 26 

showing that the displacement thickness, (5*, also decreases with increasing skin friction coefficient, c/. 

Figure 19 shows the non-dimensional boundary layer thickness, S/L, (left) and the non-dimensional dis- 
placement thickness, S* /L, (right) computed from the theoretical equations, Eqs. (24) and (26), respectively, 
for ReL = 10 TO 6 , using the local skin friction coefficient c/ of the above numerical predictions. As might have 
been deduced, the SA model yields slightly lower values of S/L and 5* / L than the SST and the SSG/LRR-w 
model, which is supposed to have an effect on the prediction of shock positions in transonic flow; see Ref. 16. 
The larger c/, the smaller the displacement thickness, d*, and the further downstream a shock is expected 
to be predicted. 


IV. A. 2. Subsonic Flow Around the NACA 0012 Airfoil 

The subsonic flow around the NACA 0012 airfoil is another basic validation test case documented on the 
TMR Website. 7 It is considered representative of mostly attached flows around airfoils and slender bodies. 
According to the TMR Website, 7 the inflow conditions are defined by a Mach number of M = 0.15 and a 
Reynolds number of Re c = 6 • 10 6 based on the chord length, c, of the airfoil. The reference temperature is 
taken as Tn e f = 300 K. 

Numerical predictions have been carried out, using CFL3D and TAU on the second finest structured 
mesh with N = 897 x 257 points provided on the TMR Website. 7 The boundary conditions employed are 
sketched in Fig. 20, also taken from the TMR Website. 7 Note that the farfield is approximately 500 chord 


9 of 47 


American Institute of Aeronautics and Astronautics 



lengths away from the airfoil, in order to minimize the effect of the lift generating surface on the farfield 
boundary condition inflow. 

According to the TMR Website, the numerical results are provided at three different incidences only, i.e., 
a = 0°, a = 10°, and a = 15°. Some complementary results have been computed with TAU at a = 1°, 
a = 2°, and a = 5°, using the SSG/LRR-w model. 

Figure 21 (left) shows the lift coefficient, Cl, versus the incidence as computed with CFL3D and TAU, 
using different turbulence models, compared to the experimental data by Ladson 22 provided on the TMR 
Website. 7 Note that all computations have been carried out in fully turbulent mode, whereas in the exper- 
iment transition was tripped. Nevertheless, according to the TMR Website, 7 these data appear to be the 
most suitable for comparison of force coefficients. 7 

Obviously, there is good agreement between all different predictions and the experiment for the lift curve, 
Fig. 21 (left). A detailed view of the results at higher incidence in Fig. 21 (right) reveals that with the SA 
model at a = 15°, the lift is predicted slightly higher than with the SST and the SSG/LRR-w model, which 
seem to be a little closer to the experimental lift curve. This difference is also observed in the force polar, 
depicted in Fig. 22, where the SA prediction might be considered in slightly better agreement with the 
measurement. However, the numerical results are too sparse for a more definite conclusion. 

The reason for the small differences between turbulence models becomes clear when looking at the 
pressure distributions at a = 10° in Fig. 23 (left) and a = 15° in Fig. 23 (right). For comparison, the 
experimental data of Gregory and O’Reilly 23 are used, which are provided on the TMR Website 7 and are 
considered most suitable, although they were obtained at a lower Reynolds number of Re c = 3 • 10 6 based on 
the chord length of the airfoil. On the overall scale there is good agreement of all numerical predictions with 
the experimental C p data, as shown in Fig. 23, whereas some differences can be observed when enlarging 
the suction peak region, as shown in Fig. 24. Clearly, around and downstream of the suction peak, the SA 
model tends to predict slightly lower pressure than the SST and the SSG/LRR-u; model. Even if there is 
still some grid influence, as can be concluded from the difference between the CFL3D and TAU results with 
the SSG/LRR-u; model, at a = 15°, this influence is obviously smaller than the difference between the SA 
and the SSG/LRR-u; predictions. 

However, in view of the experimental resolution and the uncertainty with respect to Reynolds number 
effects, any further assessment of the model predictions is impossible. All one can say is that all turbulence 
models seem to give reasonable results for this flow, where the SSG/LRR-u; and SST model predictions are 
very similar to one another and show small differences from the SA model results. 

IV. A. 3. Axisymmetric separated boundary layer 

The prediction of flow separation is a major challenge for any RANS turbulence model. As one representative 
of this class, the TMR Website 7 provides meshes and experimental data for the subsonic flow along a cylinder 
subjected to a pressure gradient by varying the contour of the surrounding outer streamline shape. The 
advantage of axi-symmetric configurations over two-dimensional configurations is that there are no possible 
side-wall effects that might spoil the flow. The flow has been studied experimentally by Driver. 24 The flow 
conditions are given by a Mach number of M = 0.08812 and a Reynolds number of ReL = 2 • 10 6 based on 
unit length of L = 1 m. The reference temperature used is Tr 6 / = 288.9 K. 

Numerical simulations have been performed using CFL3D and TAU on the second finest structured mesh 
with N = 713 x 193 points provided on the TMR Website. 7 In the simulations, the pressure gradient is 
generated by a streamline shaped outer wall with inviscid boundary conditions, where the deviation from the 
straight cylindrical wall extends from x ss —0.3m to x ss 0.85m with the maximum deflection at x ~ 0.3m. 
The grid structure and the employed boundary conditions are shown in Fig. 25. 7 

Figure 26 (left) shows the predicted pressure distributions compared to the experimental data along the 
cylinder surface in the region of the varying pressure gradient. Note that all results have been shifted in 
order to match C p = 0 at x = —0.4m. All models overestimate the pressure rise, where the SA model shows 
the largest deviation from the experimental data. Similar to the previous cases, the SST and SSG/LRR-o; 
predictions are very similar to one another. Also, TAU was run for this case with both first and second order 
turbulence advection terms. The higher order treatment makes very little difference in the results. 

The predicted distributions of the local skin friction coefficient, c/, along the surface are shown in Fig. 26 
(right). Indicated by c/ < 0, the SA model predicts a comparatively large region of separated flow, which 
is consistent with the predicted larger maximum pressure coefficient values, i.e., a larger overall adverse 
pressure gradient. The SST model predicts a slightly smaller region of separated flow, which appears to be 
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in fair agreement with the experimental data. The SSG/LRR-w model, although yielding almost the same 
pressure distribution as the SST model, predicts only a very small region of separated flow that is hardly 
visible on this scale. 

Figure 27 illustrates the flow fields predicted by the SSG/LRR-w model (left), the SA model (center) and 
the SST model (right) by depicting the corresponding streamlines in the region of the flow separation. Note 
that for clarity the wall normal y-axis is enlarged. With the SSG/LRR-ai model, the separation bubble is too 
shallow to be visible, whereas the SA and the SST models yield larger bubbles that can easily be identified. 
As will be shown in the velocity profiles, the larger bubbles are closer to the experiment. 

Figure 28 shows the corresponding profiles of the non-dimensional streamwise velocity component nor- 
malized with the inflow velocity, U c c , at six different positions in the region of the adverse pressure gradient. 
Clearly, at the most downstream position at x = —0.0762 m there is little difference between the predic- 
tions, and all models apparently agree well with the experimental data. Nevertheless, consistent with the 
observations on the pressure distribution, the SA model slightly underestimates the boundary layer edge 
velocity. In all positions further downstream, the SST model predictions maintain good agreement with the 
experiment, while in particular the SSG/LRR-w model predictions deviate from the experimental data close 
to the wall. Some very small differences between the CFL3D and the TAU results can be observed in the 
two most downstream positions, but the characteristics of the solution remain the same. 

Figure 29 shows the profiles of the non-dimensional Reynolds shear stress component, R xy , scaled with 
the square of the inflow velocity, U^, at the same six positions. At the three upstream positions, the SA 
model overpredicts the Reynolds shear stress, while it underpredicts it at the two most downstream positions. 
(Again, some differences between the SSG/LRR-w results of CFL3D and TAU can be seen at the two most 
downstream positions, but they are relatively minor.) In contrast, the SSG/LRR-cu model underpredicts the 
Reynolds shear stress at the three upstream positions and is in fair agreement with the experiment at the 
two most downstream positions. The SST model shows good agreement with the measurement at virtually 
all positions, except the last one where it underpredicts the Reynolds shear stress. 

Thus, predicted large flow separation appears to be accompanied by high upstream and low downstream 
Reynolds shear stress, and vice versa when only little separation is predicted. Overall the SST model clearly 
gives the best predictions for this case, while the small sensitivity of the SSG/LRR-w model to separation is 
a bit disappointing. 

IV. A. /. Transonic flow over an axisymmetric bump 

In transonic flow, shock induced separation may occur. As a representative of this phenomenon, the transonic 
flow along a cylinder with a bump shaped as a circular arc has been studied experimentally by Bachalo and 
Johnson. 25 The Mach number is M = 0.875 and the Reynolds number is Re c = 2.763 -10 6 based on the chord 
length, c, of the bump. The reference temperature is taken to be T^ e f = 300 K. Corresponding grids and 
experimental data are provided on the TMR Website, 7 where the grid structure and boundary conditions 
are sketched in Fig. 30. 7 The bump is located between 0 < x/c < 1, so that its crest is located at x/c = 0.5. 
Numerical simulations have been performed on the second finest structured mesh with N = 721 x 321 points. 

Figure 31 (left) shows the predicted surface pressure distributions in the range x/c = 0.4 to x/c = 1.7, 
covering the region from just upstream of the bump crest to 0.7 chord lengths downstream of the trailing 
edge of the bump. The experimental data show a suction peak at x/c ss 0.62 followed by a steep increase in 
pressure, indicating a shock. The SA model slightly overestimates the suction peak and predicts the shock 
somewhat downstream of the experimental position, whereas the SST and particularly the SSG/LRR-o; 
model predictions are in closer agreement with the experimental data. Downstream of the position x/c « 
0.75, all models tend to overpredict the pressure, where the SST model shows the smallest and the SA model 
the largest deviation. Downstream of the position x/c ss 1.4, all models predict a similar pressure coefficient, 
C p , in fair agreement with the experiment. Note that there is very little difference between the CFL3D and 
TAU predictions with the SSG/LRR-w model. 

Although no experimental data are available, the predicted skin friction distributions in terms of c/ 
versus x/c are plotted in Fig. 31 (right). Around the leading edge of the bump, the SA model predicts the 
lowest skin friction, whereas just upstream of the shock it predicts the highest c/ of all models. The SST and 
the SSG/LRR-w models yield very similar skin friction upstream of x/c « 0.1, whereas further downstream 
the SST model predicts lower c/ values that up to the shock position are the lowest of all models. In the 
experiment, shock induced separation starts at x/c « 0.7, slightly downstream of the shock position observed 
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in the pressure distribution. 25 According to Baclialo and Johnson, the flow reattaches at x/c ~ 1.1. 25 Both 
experimental separation and reattachment positions are indicated in the figure. 

Figure 31 (right) shows that all models predict separation onset directly behind the shock, indicated by 
the steep drop in c/, whereas in the experiment there is some distance in between the shock and separation. 
Consequently, the SA model predicts the position of separation onset closest to the experiment and SST the 
furthest upstream of the experimental position. Both SA and SST predict the reattachment downstream of 
the position observed in the experiment, leading to an overestimation of the length of the separation. In 
contrast, the SSG/LRR-u; model predicts both separation onset and reattachment slightly upstream of the 
respective experimental positions, leading to an overall good agreement with the experimental separation 
length. Nevertheless, the model predicts the highest c/ levels within the separated region. 

Figure 32 shows the corresponding streamline pattern as predicted by CFL3D with the different models. 
The bubbles vary considerably in length, and the SSG/LRR-u; model predicts a more rapid change in flow 
direction towards reattachment than the SA and the SST models. This behavior is probably related to the 
non-physical back-bending of streamlines observed by Jakirlic and Maduta with their model, 26 indicating 
that the SSG/LRR-u; model might benefit from a correction similar to the one described in Ref. 26. Note 
that despite the differences in c/, the SSG/LRR-u; predicts a bubble of similar thickness as the SA model, 
whereas the SST model not only yields the longest, but also the widest bubble. This observation is confirmed 
when looking at the profiles of the non-dimensional streamwise velocity, U , normalized with the freestream 
velocity, Goo, at six different streamwise positions, shown in Fig. 33. Upstream of the bump, at x/c = —0.25, 
all models predict a similar velocity profile in good agreement with the experimental data. At x/c = 0.688, 
in between the experimental shock position and separation onset, large deviations are observed between 
the different model predictions and the experiment. Close to the wall, the SA model result shows the best 
agreement with the measurement, whereas towards the boundary layer edge the SSG/LRR-u; model result 
is closest to the experiment. Nevertheless there is some deviation also with the SSG/LRR-u; model. Note 
the difference between the CFL3D and the TAU result at x/c = 0.688, which is due to high sensitivity in 
the predicted shock position. 

At positions x/c = 0.813 and x/c — 1.000, within the experimentally observed separation bubble, the 
SSG/LRR-u; model underestimates the reverse flow velocity, whereas the SST model prediction is in good 
agreement with the measurement and the SA model slightly overestimates it. Outside the bubble, the 
SA model shows the best agreement with the experimental velocity profiles, which holds for all positions 
x/c > 0.813. 

Note that at x/c = 1.000, the SSG/LRR-u; and the SA model yield zero streamwise velocity at almost 
the same position, indicating a similar thickness of the separation bubble. Apparently, the smaller absolute 
value in c/ observed with the SSG/LRR-u; is caused by a too low reverse flow speed inside the bubble, not 
because of a too small size of the separation. 

At the positions x/c = 1.125 and x/c = 1.375, downstream of the experimental reattachnrent point, the 
SSG/LRR-u; model predicts the highest velocity near the wall, which is closest to the experimental data. 
However, this might be a consequence of its prediction of a relatively small normal extent of the separation 
bubble. Note that the predicted velocity profiles obtained with the SSG/LRR-u; model and the SST model 
are very similar, except near the wall, including the separation. The SST limiter, designed to be active near 
the wall for adverse pressure gradient flows, 4 is likely the main reason for the difference. 

Figure 34 shows the corresponding profiles of the non-dimensional Reynolds shear stress, R xy , normalized 
with the squared freestream velocity, U^, at the same six streamwise positions as the above velocity profiles. 
Similar to the velocity predictions, upstream of the bump at x/c = —0.25, all model predictions are similar 
and agree well with the experimental data. At x/c = 0.688, all models show some deviations from the 
experiment, which can probably be attributed to the vicinity to the shock that is predicted at slightly different 
positions. At all positions further downstream, x/c > 0.813, the SSG/LRR-u; model always predicts the 
highest and the SA model the lowest Reynolds stresses, where the SSG/LRR-u; result is in particularly good 
agreement at x/c = 0.813 (in the separation bubble) and at x/c = 1.1375 (downstream of the reattachnrent 
point). At x/c = 1.000, which is still in the separation bubble, and at x/c = 1.125, which is in the immediate 
vicinity of the experimental reattachment point, all models underpredict the measured Reynolds shear stress. 

An assessment of models with respect to this test case appears to be particularly difficult because of two 
effects. First, the shock, as the cause of the flow separation, needs to be predicted at the correct position and 
with correct strength. Second, the model needs to react correctly to the shock, i.e., to the strong adverse 
pressure gradient in terms of the size of the separated region and accompanying details of the flow field, as 
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well as to any relaxation effects downstream of the shock. All models show certain advantages, but only 
with respect to part of the flow characteristics. 

The SA model predicts the shock downstream of the experimental position, thus the shock is stronger and 
causes a comparatively large separation bubble. Separation onset is predicted in good agreement with the 
experiment, but the bubble appears too long. This is accompanied by an underestimation of the Reynolds 
shear stress, but, except for the near-wall region, the velocity profiles agree well with the measured ones. 

The SST model apparently predicts the shock slightly upstream of the experimental position, causing 
separation onset upstream of the experimental position. The separation bubble is shorter than predicted 
by the SA model, but reattachment still occurs downstream of the experimental position. There is good 
agreement with the experimental reverse flow velocity profiles within the separation bubble, but there are 
larger deviations outside. The Reynolds shear stress is predicted larger than by the SA model, but it still 
underestimates the levels found in the experiment. 

The SSG/LRR-w model predicts the shock position in very good agreement with the experiment, but 
separation onset is found upstream of the measured position. The separation length appears to be in good 
agreement with the experiment, but this means that the reattachment point is also predicted upstream of 
the experimental position. Except for two positions, the Reynolds shear stress is predicted in rather good 
agreement with the measurements and is generally closest to the experiment compared to the other models. 
The thickness of the bubble is predicted similar to the SA model, but the reverse flow speed inside the bubble 
is underestimated. In the outer part of the boundary layer, the velocity profiles show similar deviations as 
with the SST model. 

IV. A. 5. Summary of Basic Validation Cases 

For the flow over a flat plate at zero pressure gradient and the flow around a NACA0012 airfoil up to a = 15°, 
where the flow is mainly attached, all models yield similar results in good agreement with theoretical results 
and experimental data, respectively. In the case of the separating axisymmetric boundary layer flow, the 
SSG/LRR-w underpredicts the separation extent. For the transonic flow over an axisymmetric bump, the 
SSG/LRR-w model predicts the shock location and the length of the separation bubble in (very) good 
agreement with the experiment, but the reverse flow speed inside the bubble is underestimated. 

IV. B. Extended Validation Cases 

IV.B.l. Subsonic Flow around the NACA4J12 Airfoil at a — 13.87° 

The subsonic flow around the NACA4412 airfoil at an incidence of a = 13.87° can be viewed as a simple 
case representative of high-lift flow, involving separation near the trailing edge. The flow has been studied 
experimentally by Coles and Wadcock 2 ' -28 at a comparatively low Reynolds number of Re c = 1.52 • 10 6 
based on the chord length c, where transition has been tripped at x/c = 0.025 on the upper surface and at 
x/c = 0.103 on the lower surface.' 

In contrast, according to the TMR Website, 7 the numerical simulations are carried out in fully turbulent 
mode on the second finest structured mesh with N = 897 x 257 points also provided on the TMR Website. 7 
As depicted in Fig. 35, 7 the far field is located approximately 100 chord lengths away from the airfoil. Mach 
number is M = 0.09 and reference temperature is Tfi e / = 297.8 K. 

Note that the NACA4412 experiment was in a relatively small wind tunnel, which likely had some 
influence. As a result, as discussed on the TMR website,' there are some uncertainties associated with 
running the case “correctly” vis-a-vis the experiment. Including the wind-tunnel walls might remove part of 
the uncertainties, but would considerably increase the numerical effort and probably raise questions regarding 
the prediction of secondary flow. Despite the uncertainties, this case can still be instructive for comparing 
different models’ behaviors for airfoil trailing edge separation. 

Note also that with this particular test case, a difference between CFL3D, FUN3D and TAU has been 
observed, when running the SSG/LRR-w model. While with TAU, after some initial convergence problems, 
a steady state solution has been found, with CFL3D and FUN3D the solution remains unsteady, even when 
running a time accurate simulation. With a more recent version of TAU, a different researcher also could 
not obtain a steady state solution for this case when starting from freestream conditions. However, when 
restarting from a solution at a = 10°, for which a steady state solution can be found, and successively 
increasing the incidence, the steady state solution at a = 13.87° has been confirmed. This might be an 
indication of the existence of two different solutions, one being steady and one unsteady. Therefore, only 
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steady state TAU results for SSG/LRR-w and steady CFL3D results for SA and SST are presented, and the 
steady state TAU solution needs to be considered with some caution. 

Figure 36 (left) shows the pressure distributions predicted with TAU, using the SSG/LRR-w model, 
and with CFL3D, using the SA and the SST models compared to the experimental data by Coles and 
Wadcock. 2 ' -28 The SA model predicts the highest suction peak and lower pressure than the SST and the 
SSG/LRR-w models all the way to x/c « 0.8 on the suction side, apparently in better agreement with 
the experimental data. However, the experimental data are uncorrected for wind tunnel wall effects and, 
according to the TMR Website, should be viewed in a qualitative sense, only. 7 

Figure 36 (right) shows the distribution of the local skin friction coefficient, c/, near the trailing edge 
as numerically predicted with the different models. As indicated by negative c/ values, all models predict 
separation at almost the same position, i.e., the SST model at x/c « 0.775, the SSG/LRR-w model at 
x/c ~ 0.78 and the SA model at x/c « 0.79. Interestingly, the SA and SST models yield almost identical 
values for the negative skin friction coefficient in the separation region that are almost twice as large in 
magnitude as the values predicted by the SSG/LRR-w model. 

Note that the c/ level in the reverse flow region must not be misinterpreted as a measure of the thickness 
of the separation bubble, as becomes clear from the streamline patterns depicted in Fig. 37. Indeed, the 
SST model predicts the largest separation bubble of all models tested. However, despite the different c/ 
level, the SSG/LRR-u; model predicts a separation bubble of similar thickness that encompasses a smaller 
region downstream of the trailing edge than SST. In contrast, the SA model predicts a separated region that 
is clearly shallower, although c/ is nearly identical to the SST model result. 

Figure 38 shows computed and measured profiles of the non-dimensional streamwise velocity component, 
U, along the wall normal direction at six different positions on the upper surface. Note that an unusual 
reference velocity, t/# e / « 0.93/Too, is used in the nondimensionalization of velocity and Reynolds stress. See 
the TMR website 7 for more details. Different from the results for the pressure coefficient, the SA model 
shows considerable deviation from the experimental data, predicting a too shallow separation bubble, as 
indicated by the region of negative streamwise velocity. In contrast, the predictions of the SST and the 
SSG/LRR-w models both agree quite well with the measured velocity profiles, where the SST model yields 
higher maximum reverse streamwise velocity. Note that, in particular at x/c = 0.8973 and at x/c = 0.9528, 
the SSG/LRR-w model yields the smallest and the SA model the largest absolute values of the streamwise 
velocity component in the reverse flow region, indicating that the flow in the bubble is moving at different 
speeds. This is assumed to be the major reason for the fact that the SSG/LRR-w model yields smaller 
absolute values of c/ than the SST model, although the predicted bubble sizes are apparently similar. 

Figure 39 shows the corresponding profiles of the non-dimensional Reynolds stress component R xy . All 
models underestimate the experimental values, where the SA model shows the largest and SSG/LRR-w model 
the smallest deviation from the measurement. Nonetheless, the latter difference still appears considerable. 

IV.B.2. Flow Through a Curved Duct 

The subsonic flow through a rectangular duct, undergoing a rapid 30° bend, addresses the effects of strong 
streamline curvature compared to the boundary layer thickness. The flow has been studied experimentally 
by Smits et al. 29 in a duct with a height of H = 0.127m and an aspect ratio of 6 : 1, where the focus is 
on the boundary layer on the convex (lower) wall. According to the TMR Website, 7 the flow conditions are 
defined by a subsonic Mach number of M = 0.093 and a Reynolds number of Ren = 2.667 • 10 5 based on 
the duct height H . Reference temperature is T^ef — 293.3 K. On the (upper) concave wall, considerable 
spanwise variation of the skin friction has been observed, indicating the presence of longitudinal Gortler 
vortices that may develop due to centrifugal forces. However, for the purposes of validation, the influence of 
any three-dimensionality along the upper wall is assumed to have minimal influence on the flow very near 
the convex surface, and all computations are performed in 2-D. 

The numerical predictions have been carried out on the second finest structured mesh with N = 513 x 
193 points provided on the TMR Website. 7 Figure 40 (left) shows the overall geometry together with 
the corresponding boundary conditions, while Fig. 40 (right) defines the coordinates used upstream and 
downstream of the bend. Both figures are taken from the TMR Website. 7 

Figure 41 shows the pressure distributions predicted by the different models compared to the experimental 
data. Note that all numerical results for the pressure coefficient, C p , have been shifted by A C p = —0.04, 
in order to match the experimental reference. 7 Apparently there is little difference between the model 
predictions, all agreeing well with the experimental data. 
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Larger differences are observed in the distribution of the skin friction coefficient c/ on the lower convex 
wall, Fig. 42 (left). The SSG/LRR-w result follows the admittedly sparse experimental data quite closely, 
whereas the SA model considerably overestimates the local minimum just downstream the bend. The SST 
model matches the minimum in c/, but predicts a slightly stronger increase downstream. In contrast, on the 
upper concave wall, the SSG/LRR-u; model predicts the highest skin friction of all models, as shown in Fig. 
42 (right). Nevertheless the c/ level present in the experiment is underestimated by all models. However, it 
should be noted that the simulations are carried out in two- dimensions only, preventing any Gortler vortices 
from being represented. 

Figure 43 shows the profiles of the non-dimensional wall-parallel velocity component, scaled with the 
experimental reference velocity of Un e f = 31.9m/s, at six different stations along the lower convex wall, 
where x = 0 refers to the end of the bend. As one can see, upstream of the bend at x = —0.166 ro, 
all model predictions agree well with the experimental data, whereas directly downstream of the bend at 
x = 0.030m, all models overestimate the wall parallel velocity. Further downstream, the differences become 
smaller again, where at x = 0.183m and x = 0.335?n the SSG/LRR-w model predictions appear to be in very 
good agreement with the experiment, while the SA and SST model predictions show some deviations. For 
the two most downstream positions at x = 0.635m and x = 1.250m a definite assessment is more difficult, 
since all model predictions show some deviations from the experiment. 

Concerning the corresponding Reynolds shear stress component termed R xy throughout, the SSG/LRR- 
uj model predictions are in rather good agreement with the experiment, only overestimating the minimum 
value at x = 0.355m. and x = 0.635m. At these positions, the SA and the SST models predict the minimum 
level in closer agreement with the experiment, but show larger deviations from the experiment in every other 
respect. 

Thus, there appear to be some advantages of the SSG/LRR-w model for this kind of flow. Nevertheless, 
note that SA and SST predictions can be improved for this case by employing some rotational or curvature 
correction, 7 which, however, should then be tested with other validation cases as well. 

TV. B. 3. Two-dimensional Flow Over a Wall-Mounted Hump 

The two-dimensional flow over a wall-mounted hump has been used for studying both separation and flow 
control. It has been investigated experimentally by Greenblatt et ah, 30-32 at an inflow Mach number of 
M = 0.1 and a Reynolds number of Re c = 0.936 • 10 6 based on the chord length of the hump, and served as 
a test case of the NASA CFDVAL2004 Workshop. 33 Reference temperature is Tn e f = 298.3 K. 

On the TMR Website, 7 the case is presented without any flow control in order to assess turbulence 
models’ capabilities to predict separation. Nevertheless, the geometry contains the plenum chamber inside 
the hump. The numerical predictions have been carried out on the “fine” structured mesh with N = 210, 060 
points provided on the TMR Website, 7 which is identical to one of the grids provided in the original NASA 
CFDVAL2004 Workshop. The geometry and corresponding boundary conditions are shown in Fig. 45. 7 The 
hump leading edge is located at x/c= 0. 

Note that for this particular test case, only results with CFL3D and FUN3D are presented, as considerable 
differences have been observed between the TAU and the CFL3D and FUN3D SSG/LRR-w results to date. 
Additional tests with TAU revealed a sensitivity of the solution to a variation in the outlet boundary 
condition; however, the reason for the differences has not yet been identified. Therefore the TAU results 
have been omitted. 

Figure 46 (left) shows the computed distributions of the local pressure coefficient, C p , along the wall 
compared to the experimental data by Greenblatt et al. 30-32 Note that all simulation data have been shifted 
by A C p = —0.015, in order to better match the experimental reference upstream. All models underpredict 
the suction peak at x/c « 0.5. Further downstream, the experimental data show a local maximum in C p 
at x/c k. 0.7, followed by another small suction peak at x/c « 0.85. Apparently, none of the models can 
represent this behavior in close agreement with the measurement. The SA model gives the correct shape, 
but overpredicts the local pressure maximum at x/c « 0.7. In contrast, the SST and the SSG/LRR-w 
models predict the local pressure maximum in good agreement with the experiment, but in particular the 
SSG/LRR-w fails to predict the second suction peak. Downstream of the second suction peak, the SA model 
is in fair agreement with the experimental data, whereas the SST and the SSG/LRR-w models show larger 
deviations. Note that there is no noticeable difference between the CFL3D and the FUN3D results with the 
SSG/LRR-w model. 
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Obviously, the behavior in the range 0.7 < x/c < 1.2 is related to flow separation, as indicated in Fig. 46 
by the region of negative skin friction coefficient, c/. The SSG/LRR-w model predicts the streamwise extent 
of the separation bubble in fair agreement with the experiment and consequently shows the least deviations 
in the downstream recovery region. However, it underestimates the c/ levels in the separation bubble that 
are predicted by the SST and the SA models in much better agreement with the experiment. In contrast, 
the SA and SST models overestimate the streamwise extent of the separation and show larger deviations in 
the recovery region. 

Figure 47 illustrates the predicted flow fields in terms of the respective streamline patterns. As indicated 
by the skin friction distribution, the SSG/LRR-ta model predicts the shortest bubble, but apparently also 
the widest one. As with the transonic axisymmetric bump case, the SSG/LRR-o; model predicts bubble 
reattachment with a more rapid change in flow direction (i.e., “back-bending,” as discussed in section IV. A. 4). 
The SA and the SST models predict a longer and, in particular the SA model, a shallower bubble. Apparently, 
the low absolute value of the skin friction coefficient in the separation region predicted by the SSG/LRR-o; 
model is an indication of too low reverse flow speed, but is not an indication of too little separation. Thus, 
in this case, the apparently improved agreement with the experimental Cf data near to and downstream of 
the reattachment point should probably not be interpreted as a compensation for insensitivity of the model 
to separation. 

This is confirmed when considering the velocity profiles. Figure 48 shows profiles of the non-dimensional 
streamwise velocity component, U, scaled with the reference velocity, U^ef, at six different stations along 
the bubble, beginning at x/c = 0.80, where the flow is already separated. Clearly, at the three positions 
x/c = 0.80, x/c = 0.90 and x/c — 1.00, the SSG/LRR-w model underestimates the reverse flow velocity, and 
shows the widest region of reverse flow. Apparently, it is the too-small reverse flow speed that is responsible 
for the too low absolute values of c/ observed in Fig. 46. The SST model shows the best agreement with the 
measured data of the reverse streamwise flow component, but also considerably overestimates the width of 
the reverse flow region. The SA model tends to overestimate the maximum reverse flow speed, but at these 
three positions is in very good overall agreement with the experimental data. In this respect, the situation 
resembles the observation with the separated axisymmetric boundary layer in section IV. A. 3. One has to 
conclude that the experimental separation bubble is indeed shallower than predicted by the SSG/LRR-w 
and the SST models. Note that the CFL3D and FUN3D results with the SSG/LRR-w model predict nearly 
identical velocity profiles. 

Figure 49 shows the corresponding profiles of the non-dimensional Reynolds stress component, R xy , 
scaled with the square of the reference velocity, U^ e j. All models tend to underestimate the level of the 
Reynolds stress R xy , where, except at position x/c = 0.80, the SSG/LRR-w gives the highest negative values. 
Apparently, higher upstream Reynolds stress levels are associated with a shorter separation bubble, but there 
is no obvious relation to its width. 

In the two most downstream positions, the SSG/LRR-w model predicts the profiles of the non-dimensional 
Reynolds stress component, R xy /U^p in fair agreement with the experimental data. This is probably a 
consequence of the more upstream reattachnrent point, which, according to the distribution of the skin 
friction coefficient c/ shown in Fig. 46 (right), is in closer agreement with the experiment than with the 
other models. 

The differences between the CFL3D and FUN3D results with the SSG/LRR-w are somewhat larger for the 
non-dimensional Reynolds stress component, R xy /D\ e f, than for the other quantities shown. Nevertheless 
they are clearly smaller than the differences between the results for different models. 

Thus, all models show a combination of advantages and disadvantages. The SSG/LRR-w model predicts 
a separation bubble that has approximately the correct length, but is too thick, and the reverse flow speed 
inside the bubble is too low. The SST model predicts a separation bubble that is too wide and too long, 
whereas the reverse flow speed agrees with the experiment. The SA model predicts a separation bubble that 
probably has the correct width, but is too long and the reverse flow speed is higher than in the experiment. 
Furthermore, there appears to be only a relation between the predicted level of the Reynolds stress component 
R xy and the predicted length of the separation bubble, but not with its width. 

IV. B./. Supersonic flow Through a Square Duct 

In a duct with rectangular cross section, secondary flows develop longitudinal vortices that are caused by 
the difference between the Reynolds normal stresses close to walls. This phenomenon cannot be represented 
by linear eddy-viscosity models due to the underlying Boussinesq hypothesis. 
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The corresponding test case on the TMR Website' involves the supersonic flow at a Mach number of 
M = 3.9 through a square duct of height and width D = 25.7 mm over a length of L = 50 D. The Reynolds 
number is Re d = 5.08 • 10 5 based on the duct height D. Reference temperature is taken to be = 288.9 
K. The flow has been studied experimentally by Davis and Gessner. 34 Figure 51 (left) shows an overall 
sketch of the geometry together with the corresponding boundary conditions. Figure 51 (right) shows the 
grid structure in a cross-section of the duct. Both figures are taken from the TMR Website.' Only one 
quarter of the duct is simulated. 

The numerical predictions are typically carried out on the second finest structured mesh provided on the 
TMR Website 7 with N = 481 x 81 x 81 points. However, for this case it was found that the deviations 
between the predictions of CFL3D and TAU on the second finest grid using the SSG/LRR-w are larger than 
observed with the other test cases shown previously. For this reason, FUN3D results have been included, 
confirming the larger spread. These differences indicate a larger grid influence with the SSG/LRR.-w than 
with the eddy-viscosity models for this case. This is also confirmed by computations with TAU on the finest 
grid with N = 961 x 161 x 161 points. 

Figure 51 shows the distribution of the local skin friction coefficient along the vertical wall at xj D = 50 
as predicted by the different models compared to the experimental data. Most of the CFD results are from 
the second-finest grid, but the TAU results are also plotted for SSG/LRR-w on the finest grid. There is 
some small change due to grid refinement. The SSG/LRR-w model yields the skin friction coefficient in fair 
agreement with the experimental data, except for the duct center. In contrast, the SA and SST models 
results show large deviations throughout, which can be attributed to the missing secondary flow they cannot 
predict. 

Figures 52 and 53 show the numerical results for the non-dimensional streamwise velocity component 
along the vertical centerline (left) and along the diagonal (right) at xj H = 40 and 50 compared to the respec- 
tive experimental data. The SSG/LRR-ai model predicts profiles in fair agreement with the measurement. 
In contrast, the SA and the SST model tend to overpredict the non-dimensional streamwise velocity along 
the vertical centerline, whereas they underpredict it along the diagonal. As already stated, this is expected 
due to the general inability of linear eddy-viscosity models to predict the secondary flow. 

Figure 54 shows the velocity contours in the cross-section at x/H = 40 as predicted by the three codes, 
using the SSG/LRR-u; model. The two counter-rotating longitudinal vortices are clearly visible. Although 
not shown here, the linear eddy-viscosity models show no indication of these structures. For comparison, 
the TAU results on the finest grid are also depicted, illustrating that the results for the SSG/LRR-oi model 
on the second finest mesh are not yet grid converged. 

IV. B. 5. Summary of Extended Validation Cases 

For the flow around the NACA4412 airfoil, the SSG/LRR-u; model predicts the trailing-edge separation 
similar to the SST model, but yields a lower reverse flow velocity. For the flow over a wall-mounted hump, the 
SSG/LRR.-w model predicts the length of the separation bubble in fair agreement with the experimental data, 
but apparently overestimates its width. The reverse flow velocity inside the bubble is again underestimated. 
For the curved duct case the SSG/LRR-u; model yields the best results of all models, where in particular the 
velocity profiles downstream of the bend are in very good agreement with the experiment. For the supersonic 
flow through a square duct, the SSG/LRR-u; model predicts secondary flow and shows fair agreement with 
the experimental skin friction and velocity data in two downstream cross-sections. As expected, for the latter 
two flow cases, the SSG/LRR-u; model is superior to the SA and the SST models. 

V. Conclusion 

The SSG/LRR-u; differential Reynolds stress model has been implemented into the flow solvers CFL3D 
and FUN3D by NASA and into TAU by DLR. The respective implementations have been cross-checked 
by running three verification test cases of the TMR Website 7 on a series of successively refined structured 
grids: the flow over a flat plate with zero pressure gradient, a plane jet flow and the two-dimensional 
flow over a hump. The development of characteristic integral and local quantities with increasing mesh 
resolution has been checked, showing that the different results apparently converge to a common solution. 
Excellent agreement is found, comparing the pressure and the skin friction distribution as well as velocity 
and Reynolds stress profiles obtained on the respective finest grid, giving confidence in the correctness of 
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the model implementation into the different codes. 

A series of four basic validation test cases has been studied running CFL3D and TAU with the SSG/LRR- 
oj model: the flow over a flat plate at zero pressure gradient, the flow around the NACA0012 airfoil at different 
angles of attack, an axisymmetric separating boundary layer flow and the transonic flow over an axisymmetric 
bump. The results have been compared against experimental data as well as simulation data obtained with 
CFL3D, using the SA and the SST models. 

For the flat plate flow, all results are very similar and in good agreement with experimental skin friction 
data and the log-law of the velocity profile. It is noted that the SSG/LRR-u; and the SST model results are 
almost identical. The SA model predicts slightly higher skin friction, which is probably associated with a 
slightly smaller boundary layer thickness, 5, and displacement thickness 5* at the same location. 

Good agreement between the models is also found with the NACA0012 case, where, only at a = 15°, the 
SSG/LRR-w and the SST models yield slightly lower lift and higher drag than the SA model. Apparently 
this is related to a higher suction peak predicted by the SA model. 

Similarly, for the separated axisymmetric boundary layer case the SSG/LRR-w and the SST models 
predict a smaller pressure rise in better agreement with the experiment than the SA model. Nevertheless, the 
SSG/LRR-w model significantly underestimates the associated separation, whereas the SST model predicts 
the flow characteristics in good agreement with the experimental data. With the SA model, the length of 
the separation is slightly overpredicted, but the velocity profiles are in good agreement with the measured 
ones. 

The transonic flow over an axisymmetric bump shows the largest differences between the models and is 
the basic test case that is the most difficult to be assessed. The reason is that the models predict the shock 
at slightly different positions and thus, with different strength. The SSG/LRR-w model predicts the shock 
position and the length of the separation bubble in good agreement with the experiment, but the reverse flow 
speed inside the bubble is underestimated. This is probably the cause for too low absolute values of the skin 
friction coefficient in the reverse flow region. The SST model predicts the shock slightly upstream, whereas 
the SA model predicts it downstream of the experimental position. Both models yield longer separation 
bubbles and larger absolute values of the skin friction coefficient in the reverse flow region. Downstream of 
the shock the SA model shows the overall best agreement with the experimental velocity profiles. 

Finally, four different flows from the list of extended validation test cases on the TMR Website 7 have 
been studied: separated flow over a NACA4412 airfoil, the flow through a curved duct, separated flow over 
a wall-mounted hump, and supersonic flow in a square duct. 

The flow around the NACA4412 airfoil at incidence a = 13.87° involves a trailing-edge separation. This 
case is difficult to compute with SSG/LRR-w, and there is some indication of the existence of multiple 
solutions, one steady and one unsteady. The steady SSG/LRR-w solution and the SST solution yield 
similarly large bubbles, where the velocity profiles indicate fair agreement with the experiment. In contrast, 
the bubble predicted by the SA model is smaller. Once again the smaller absolute value of the skin friction 
coefficient predicted by the SSG/LRR-w model compared to the SST model appears to be due to a lower 
reverse flow speed in the bubble. 

The flow through a curved duct addresses effects of large streamline curvature that cannot be expected 
to be represented well by any eddy-viscosity model. In contrast, the SSG/LRR-w predicts the pressure and 
skin friction distributions and even profiles of the Reynolds shear stress in rather good agreement with the 
experiment. 

The flow over a wall mounted hump again involves a large region of separated flow that is predicted by 
all models differently. The SSG/LRR-w models predict the length of the bubble in fair agreement with the 
experiment, but its width appears too large. With the SA model, the bubble width is probably close to 
experiment, but reattachnrent is too far downstream. Finally, the SST model apparently overestimates both 
length and width of the separated region. The SSG/LRR-w model yields the best agreement with respect 
to Reynolds stress profiles. 

The supersonic flow through a square duct is characterized by secondary flow that cannot be represented 
by linear eddy-viscosity models. In contrast, the SSG/LRR-w model predicts the associated longitudinal 
vortices and represents skin friction and velocity profiles in fair agreement with the experimental data. 

To summarize: for a wide range of separated flows the SSG/LRR-w model yields similar results as 
the SA or the SST models. Advantages have been observed with respect to the prediction of the shock 
position and the length of separation bubbles. However, the reverse flow speed inside the bubble seems to 
be generally underestimated. Furthermore, effects of streamline curvature and secondary flow are predicted 
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reasonably well by the SSG/LRR-ct; model. Capturing these effects requires curvature corrections or non- 
linear extensions in standard eddy- viscosity models. Thus, the SSG/LRR-o; model appears to have a wider 
range of applicability than eddy-viscosity models without special modifications. 

Future work on the SSG/LRR-u; model will be needed in order to improve the prediction of reverse flow 
speed in separation bubbles, and to eliminate the “back-bending” that occurs near bubble reattachment. 
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Figure 1. Flow over a flat plate at zero pressure gradient. Boundary conditions according to TMR Website. 7 
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Figure 2. Flow over a flat plate at zero pressure gradient. Grid convergence of total drag coefficient (left) and local 
skin friction coefficient at x/L = 0.485 (right). 2nd order refers to turbulence advection treatment. 



Figure 3. Flow over a flat plate at zero pressure gradient. Predicted skin friction distributions on the finest mesh with 
N = 545 X 385 grid points. 


21 of 47 


American Institute of Aeronautics and Astronautics 



Figure 4. Flow over a flat plate at zero pressure gradient. Predicted velocity profiles in inner scaling at position 
xjh = 0.485. Results for finest mesh with N = 545 X 385 grid points. 




Figure 5. Flow over a flat plate at zero pressure gradient. Predicted Reynolds stress profiles in inner scaling at position 
x/L — 0.485. Results for finest mesh with N = 545 X 385 grid points. 


22 of 47 


American Institute of Aeronautics and Astronautics 


1 20 r- 

100 : 
80 : 

_ 60 - 

>. 

40 - 

20 : 
0 : 
-20 : 


symmetry 


Pt/P rpf =1 .046 
-Tt/T re( =1 .0 

1 quantity from interior 


P/P re| =1.0, 

other quantities 
from interior 


y see close-up plot 


50 


symmetry 


100 

x/H 


150 


200 


1.5 

1 

C 

>> 

0.5 

0 

-0.5 


Pt/P„, =1.046 
Tt/T [el =1 .0 

1 quantity from interior 


adiabatic solid wall 


/ 


Pt/P .=1.1862 
Tt/T rel =1 .05 

1 quantity from interior 


V 

adiabatic solid wall 


_ | Close-up] 

_J I !_ 

-10 


symmetry 


-5 


x/H 


Figure 6. Plane jet developing behind an orifice. Boundary conditions according to TMR Website. 7 




Figure 7. Plane jet developing behind an orifice. Grid convergence of total drag coefficient on the plate (left) and the 
non-dimensional streamwise velocity velocity Uo/aR e f on the jet centerline (right). 



Figure 8. Plane jet developing behind an orifice. Non-dimensional streamwise velocity Uo/aR e f on the jet centerline 
computed on the finest mesh with N = 327, 680 points. 
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Figure 9. Plane jet developing behind an orifice. Profile of non-dimensional streamwise velocity across the jet at 
position x/H — 29.2468. Results on finest mesh with N — 327,680 points. 




Figure 10. Plane jet developing behind an orifice. Profiles of non-dimensional Reynolds stresses across the jet at 
position x/H — 29.2468 (R a a represents R X xi Ryyi or Rzz)- Results on finest mesh with N — 327,680 points. 
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Figure 11. Two-dimensional bump in a channel. Boundary conditions according to TMR ’Website. 7 




Figure 12. Two-dimensional bump in a channel. Grid convergence of lift coefficient (left) and viscous and pressure 
drag coefficient (right). 
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Figure 13. Two-dimensional bump in a channel. Grid convergence of local skin friction coefficient at three different 
stations. 
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Figure 14. Two-dimensional bump in a channel. Pressure distribution (left) and skin friction distribution (right) 
computed on the finest mesh with N = 1,409 X 641 points. 



Figure 15. Two-dimensional bump in a channel. Velocity profile in inner scaling at position x/H — 24. Results on the 
finest mesh with N — 1,409 X 641 points. 




Figure 16. Two-dimensional bump in a channel. Reynolds stress profiles in inner scaling at position x/H = 24 (R 
represents R xx , or Results on the finest mesh with N = 1,409 X 641 points. 
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Figure 17. Flow over a flat plate at zero pressure gradient. Velocity profile in inner scaling. Comparison of numerical 
predictions on the finest mesh with N = 545 X 385 points with analytical descriptions of the viscous sublayer and the 
log- law. 




Figure 18. Flow over a flat plate at zero pressure gradient. Skin friction coefficient vs. local Reynolds number. 
Comparison of numerical results on the finest mesh with N = 545 X 385 points with experiments by Wieghardt 18 as 
provided by Coles and Hirst. 20 




Figure 19. Flow over a flat plate at zero pressure gradient. Non-dimensional boundary layer thickness according to 
theoretical equation (24) (left) and displacement thickness according to theoretical equation (26) (right) for Rcl — 10- 10 6 
using numerically predicted ^-distributions. 
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Figure 20. Subsonic flow around NACA 0012 airfoil. Boundary conditions according to TMR Website. 7 
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Figure 21. Subsonic flow around NACA 0012 airfoil. Lift coefficient versus incidence. Comparison of numerical results 
on the second finest grid with N — 897 X 257 points with experiments by Ladson. 22 



Figure 22. Subsonic flow around NACA 0012 airfoil. Force polar. Comparison of numerical results on the second finest 
grid with N = 897 X 257 points with experiments by Ladson. 22 
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Figure 23 . Subsonic flow around NACA 0012 airfoil. Pressure distribution at a = 10° and ct — 15°. Comparison of 
numerical results on the second finest grid with N — 897 X 257 points with experiments by Gregory and O’Reilly. 23 




Figure 24 . Subsonic flow around NACA 0012 airfoil. Detail of the pressure distribution at ot = 10° and ol = 15° near 
the suction peak. Comparison of numerical results on the second finest grid with N — 897 X 257 points with experiments 
by Gregory and O’Reilly. 23 
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Figure 25. Axisymmetric separated boundary layer. Grid structure and boundary conditions according to TMR 
Website. 7 Dash-dot ellipses represent the cylinder. 




Figure 26. Axisymmetric separated boundary layer. Pressure distribution (left) and skin friction distribution (right). 
Comparison of numerical results on the second finest mesh with N = 713 X 193 points with experiments by Driver. 24 
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Figure 27. Axisymmetric separated boundary layer. Streamlines in the region of flow separation predicted with CFL3D 
using the SSG/LRR-u; model (left), the SA model (center) and the SST model (right). 
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Figure 28. Axisymmetric separated boundary layer. Non-dimensional velocity profiles at different streamwise positions. 
y is the wall normal coordinate originating from the surface. Comparison of numerical results on the second finest 
mesh with N = 713 X 193 points with experiments. 24 
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Figure 29. Axisymmetric separated boundary layer. Profiles of non-dimensional Reynolds shear stress at different 
streamwise positions, y is the distance from the surface. Reynolds stresses multiplied by a factor of 1000. Comparison 
of numerical results on the second finest mesh with N — 713 X 193 points with experiments. 24 
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Figure 30. Transonic flow over an axisymmetric bump. Grid structure and boundary conditions according to TMR 
Website. 7 Dash-dot ellipses represent the cylinder. 




Figure 31. Transonic flow over an axisymmetric bump. Pressure distribution (left) and skin friction distribution (right). 
Comparison of numerical results on the second finest mesh with N = 721 X 321 points with C p data from experiments of 
Bachalo and Johnson. 25 Experimental separation and reattachment indicated in Cf plot by small vertical lines. 



Figure 32. Transonic flow over an axisymmetric bump. Streamlines in the region of flow separation predicted with 
CFL3D using the SSG/LRR-cj model (left), the SA model (center) and the SST model (right). 
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Figure 33. Transonic flow over an axisymmetric bump. Non-dimensional velocity profiles at different streamwise 
positions, y is the coordinate normal to the streamwise direction, originating from the surface. Comparison of 
numerical results on the second finest mesh with N = 721 X 321 points with experiments by Bachalo and Johnson. 25 
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Figure 34. Transonic flow over an axisymmetric bump. Non-dimensional Reynolds shear stress profiles at different 
streamwise positions, y is the coordinate normal to the streamwise direction, originating from the surface. Comparison 
of numerical results on the second finest mesh with N — 721 X 321 points with experiments by Bachalo and Johnson. 25 
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Figure 35. Subsonic flow around NACA 4412 airfoil. Boundary conditions and coordinate system, according to TMR 
Website. 7 




Figure 36. Subsonic flow around NACA 4412 airfoil. Pressure distribution (left) and skin friction distribution near 
trailing edge (right). Comparison of numerical results on the second finest mesh with N — 897 X 257 points with C p data 
from experiments of Coles and Wadcock. 27, 28 

I SSG/LRR-oo, TAUl ISA, CFL3D1 |SST, CFL3D| 





Figure 37. Subsonic flow around NACA 4412 airfoil. Streamlines in the region of flow separation predicted with TAU 
using the SSG/LRR-u; model (left) and with CFL3D, using the SA model (center) and the SST model (right). 
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Figure 38. Subsonic flow around NACA 4412 airfoil. Non-dimensional streamwise velocity profiles at different stream- 
wise positions, d is the coordinate normal to the wall originating from the upper surface. Comparison of numerical 
results on the second finest mesh with N = 897 X 257 points with experiments by Coles and Wadcock. 27 ’ 28 
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Figure 39. Subsonic flow around NACA 4412 airfoil. Profiles of non-dimensional Reynolds stress component R xy at 
different streamwise positions, d is the coordinate normal to the wall originating from the upper surface. Comparison 
of numerical results on the second finest mesh with N = 897 X 257 points with experiments by Coles and Wadcock. 27 ’ 28 
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Figure 40. 


Flow through a curved duct. 


Boundary conditions and coordinate system, according to TMR Website. 7 



x in [m] 


Figure 41. Flow through a curved duct. Pressure distribution in terms of pressure coefficient, C p , on bottom wall. 
Comparison of numerical results on the second finest mesh with N = 513 X 193 points with experiments by Smits et al. 29 
According to TMR Website, 7 the numerical results are shifted by A C p = —0.04, in order to match the experimental 
reference. 

| Bottom (convex) waif) |Top (concave) wall | 




Figure 42. Flow through a curved duct. Skin friction distribution in terms of local skin friction coefficient on bottom 
wall (left) and top wall (right). Comparison of numerical results on the second finest mesh with N = 513 X 193 points 
with experiments by Smits et al. 29 
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Figure 43. Flow through a curved duct. Non-dimensional velocity profiles at different streamwise positions, d is the 
coordinate normal to the wall originating from the lower surface (see Fig. 40 (right)). Comparison of numerical results 
on the second finest mesh with N = 513 X 193 points with experiments by Smits et al. 29 
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Figure 44. Flow through a curved duct. Non-dimensional Reynolds shear stress profiles at different streamwise 
positions, d is the coordinate normal to the wall originating from the lower surface (see Fig. 40 (right)). Comparison 
of numerical results on the second finest mesh with N — 513 X 193 points with experiments by Smits et al. 29 
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Figure 45. Two-dimensional flow over a wall-mounted hump. Boundary conditions according to TMR Website. 7 




Figure 46. Two-dimensional flow over a wall-mounted hump. Pressure distribution (left) and skin friction distribution 
(right). Comparison of simulation results on the fine mesh with N — 210,060 points with experiments by Greenblatt et 
al. 30 " 32 
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Figure 47. Two-dimensional flow over a wall-mounted hump. Streamlines in the region of flow separation predicted 
with CFL3D using the SSG/LFtR-u; model (left), the SA model (center) and the SST model (right). 
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Figure 48. Two-dimensional flow over a wall-mounted hump. Profiles of the non-dimensional streamwise velocity 
component at different positions long the hump. Comparison of numerical results on the fine mesh with N = 210, 060 
points with experiments by Greenblatt et al. 30-32 
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Figure 49. Two-dimensional flow over a wall-mounted hump. Profiles of the non-dimensional Reynolds stress component 
R xy at different positions long the hump. Comparison of numerical results on the fine mesh with N = 210, 060 points 
with experiments by Greenblatt et al. 30-32 
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Figure 50. Supersonic flow through a square duct. Boundary conditions and grid structure, according to TMR Website. 7 
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Figure 51. Supersonic flow through a square duct. Skin friction distribution on vertical wall. Comparison of numerical 
results on the second finest mesh with N = 481 X 81 X 81 points with experiments by Davis and Gessner. 34 “SSG/LRR-o;, 
TAU, fine” refers to TAU result, using SSG/LRR-w on the finest mesh, with N = 961 X 161 X 161 points. 
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Figure 52. Supersonic flow through a square duct. Non-dimensional velocity at position x/D = 40 along vertical 
centerline (left) and along the diagonal (right). Comparison of numerical results on the second finest mesh with 
N = 481 X 81 X 81 points with experiments by Davis and Gessner. 34 “SSG/LRR-w, TAU, fine” refers to TAU result, 
using SSG/LRR-co on the finest mesh with N — 961 X 161 X 161 points. 
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Figure 53. Supersonic flow through a square duct. Non-dimensional velocity at position x/D — 50 along vertical 
centerline (left) and along the diagonal (right). Comparison of numerical results on the second finest mesh with 
N — 481 X 81 X 81 points with experiments by Davis and Gessner. 34 “SSG/LRR-o), TAU, fine” refers to TAU result, 
using SSG/LRR-o; on the finest mesh with N = 961 X 161 X 161 points. 
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Figure 54. Supersonic flow through a square duct. Velocity field at x/H = 40. Numerical results with CFL3D, 
FUN3D and TAU on the second finest mesh with N = 481 X 81 X 81 points and with TAU on the finest mesh with 
N — 961 X 161 X 161 points, using the SSG/LRR-w model. Colors refer to absolute velocity in the plane, scaled with the 
streamwise centerline velocity. 
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